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Abstract: We present the QCD analysis of /^-parameter distribution in 
near-to-planar 3-jet e^e~ annihilation events. We derive the all-order resummed 
perturbative prediction and the leading power suppressed non-perturbative correc- 
tions both to the mean value and the distribution. Here non-perturbative corrections 
are larger than in 2-jet shape observables, so that higher order non-perturbative ef- 
fects could be relevant. Experimental data (not yet available) are needed in order 
to cast light on this important point. The technique we develop aims at improving 
the accuracy of the theoretical description of multi-jet ensembles, in particular in 
hadron-hadron collisions. 
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1. Introduction 

The analysis of event shape distributions in e+e^ has provided various tests of QCD 
[1] and measurements of the running coupling [2] . The shape observables which have 
been most intensively studied and tested are: thrust T, C-parameter, jet mass 
and broadening B. Their distributions are collinear and infrared safe (CIS) and 
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therefore can be computed order-by-order in perturbation (PT) theory. Specially 
interesting is the region of two narrow jets {1 —T,C, M'^ /Q"^, B ^ 1), where the 
PT expansion needs to be resummed so that the QCD structure is most intensively 
probed. The available PT results [3] for these typical 2-jet observables^ involve all- 
order resummation of double- (DL) and single-logarithmic (SL) contributions and 
matching of the approximate resummed expressions with the exact second order 
matrix elements. 

To make quantitative predictions one needs to go beyond PT calculations and to 
take into account the l/Q-suppressed power corrections arising from the interaction 
in the confinement region. It has been proposed [4] that these non-perturbative 
(NP) corrections can be estimated by extrapolating the running coupling into the 
large distance region. A systematic method for this extrapolation is provided by 
the dispersive approach [5]. The leading NP corrections to the shape distributions 
involve a single parameter, usually denoted by ao, which is given by the integral of 
the QCD coupling over the region of small momenta k < /ij (with the infrared scale 
fjLj conventionally chosen to he ni — 2 GeV) . They have been computed at two- loop 
level to take into account effects of non-inclusiveness of jet observables [6]. Such 
a procedure for computing the l/Q-power corrections is consistent with the data 
and the NP parameter ao has been measured and appears to be universal with a 
reasonable accuracy [1]. 

We have recently studied the distribution in the thrust minor T„j, a CIS shape 
observable characterising 3-jet events. It starts at order and measures the radia- 
tion out of the event plane in e^e~ annihilation. The event plane is defined by the 
T and Tm axes, with Tm is the thrust major. The most interesting region in which 
one probes QCD is that of nearly planar 3-jet events, <^ Tm ~ T, where PT 
resummation is needed. The methods developed for the analysis of 2-jet observables 
T, C, M^/Q"^ and B were extended to the case of Tm in [7]. 

In this paper we report the analysis of the distribution in another typical 3- 
jet observable, the D-parameter which has been introduced in [8]. This is a CIS 
shape observable which, as T„i, measures the radiation out of the event plane. In 
the nearly planar 3-jet region {D <^ 1), the main difference between D and 
is that only soft particles at large angles with respect to the event plane, pout ~ 
Pin o: D <^ 1, contribute to D, contrary to the case. The fact that energetic 
hadrons with p^^ ^ pout do contribute to gives rise, in particular, to additional 
logarithmic enhancement of the NP contribution from the regions where a small 
transverse momentum energetic gluon is emitted close to the direction of one of the 
three jets. This difference is similar to that between B and T,C,M'^ /Q^ in 2-jet 

^Sometime these observables are called 3-jet observables since the first contribution involves 
three particles in the final state. We prefer to call them 2-jet observables since we are interested to 
the kinematical region in which one particle of the three is soft. 
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events [6]: an energetic hadron emitted near the thrust axis contributes to B but 
does not contribute to the other shape observables. 

Our study is performed with the accuracy necessary to make quantitative pre- 
dictions. We perform the all-order resummation of DL- and SL-enhanced PT con- 
tributions, match it [3] with the exact fixed order matrix element calculation and 
compute the leading 1/Q NP correction at the two-loop order. Actually, the fact 
that only hadrons at large angles contribute to D <^ 1 makes the present analysis 
significantly simpler than that of T^. 

We find that, as in the case, the structure of the result is quite rich, especially 
for small D where both the PT and NP components of the D distribution essentially 
depend on the geometry of the event (the angles between jets). Not only will the 
analysis of the D distribution provide an alternative measurement of the QCD cou- 
pling ctg. It gives a powerful tool for accessing genuine confinement hadronisation 
effects, for extracting the NP parameter aQ and testing its universality. 

The paper is organised as follows. In section 2 we introduce the observable and 
discuss the characteristic kinematics in the near-to-planar 3-jet region. In section 3 
we discuss the resummation procedure. In section 4 we analyse the PT resummation 
at SL level. In section 5 we analyse the leading NP corrections in terms of the uni- 
versal parameter oq. In section 6 we report the final result and discuss the necessary 
ingredients of the one-loop matching. In section 7 we report the numerical analysis 
of the distribution and the mean. In section 8 we summarise and discuss the results. 



2. The observable 

The D parameter is defined as [8] 

D ^27 det 6 = 27X^X2X3, ^a/3 = ^^4^ET^' (2-1) 

T,h\Ph\^ \Ph\ 

with the momentum of the emitted hadron h and Pha its a-component (a = 1 , 2, 3) . 
The eigenvalues A, satisfy the conditions < Aj < 1 and Aj = 1. We assume the 
order A3 < A2 < Ai. 

To select near-to-planar events we introduce a lower limit yc for the 3-jet reso- 
lution variable ys, defined^ according to the kx (Durham) algorithm [10]. We study 

^Given the set of all momenta, one defines the "jettiness" variable [9] of any pair ph and ph' 
as the quantity yhh' = 2 (1 — cosOhh') iniii(i^^, The pair of momenta pj^Phi with the 
minimum distance 2/3 — i/h^j^i = min/j^/ yhh' are substituted with the pseudoparticle (jet) momentum 
Ph" — Ph ~^Ph'- procedure is repeated with the new momentum set till only three jets are left. 
Then the final value of 2/3 is defined as the three-jet resolution of the event. 
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the following normalised integrated distribution 



E(L', y,) = V / dam 6(1/3 - Vc) 6 (L'-27det 9) , 



(2.2) 



where denotes the differential distribution in the m final particles with mo- 
menta ph- 

We comment here on the fact that we use the variable to select the 3-jet 
events (instead of, for instance, the thrust T), which makes it easier to interpret 
the 1/Q power corrections that, as we shall see, are present in the D-distribution. 
Indeed, the difference between the hadron and parton level values of the thrust 
variable T is known to be of the order of while the values of ^3 for hadrons and 
partons is not affected by l/Q corrections [5]. Therefore, the 1/Q contribution to 
the distribution (2.2) can be looked upon as a genuine NP correction to the variable 
D. On the contrary, substituting T < Tc for y^ > yc in the distribution (2.2), the 
1/Q correction would have to be shared between D and T. 

The region of small D and relatively large yc (we consider typical values of yc in 
the range 0.025 — 0.1) corresponds to the region 

A3 < A2 < Ai , Tm<^TM<T. (2.3) 

Here Ai,A2,y3 and T,Tm are determined by the hard momenta characterising 
the three jets. The smallest eigenvalue A3 and are determined by the out-of- 
event-plane soft momentum component of the particles around and between the jets 
(inter- and intra-jet radiation). Taking the plane of the event as the {y, 2;}-plane, 
the dominant contribution to A3 in the region (2.3) is given by 

pI 



Aa-EfTi- (2.4) 



EhQ 

Due to the energy factor Eh in the denominator, A3 gets the leading contribution 
from hadrons with Phx- Instead, the observable T^, 

= E ^ ' (2-5) 
h ^ 

as mentioned in the Introduction, receives contributions from particles with arbitrary 
large energies, and in particular those that are quasi-coUinear with one of the hard 
jets. 
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3. Part on process and resummation 



Our first aim is to express the distribution 'E{D,yc) in terms of parton processes in 
tlie near-to-pianar region (2.3). Here tlie parton events can be treated as 3-jet events 
generated by a hard quark- antiquark-gluon system accompanied by an ensemble 
of secondary partons ki. In this region, the quantities Ai, A2, T, Tm and the 3- 
jet resolution variable j/3 are determined by the three momenta of the hard quark, 
antiquark and gluon which we denote by Pi, P2 and P3, respectively. Introducing the 
Born x-invariants (for the Born Kinematics see Appendix A) 

^"=~^^' 2^Xi+X2+X3, (3.1) 

we have 

2(1 -a;i)(l -a;2)(l -xs) Xmm{l - x^a^) , . 

M M — , — ~ ' ^ ' 

X-1X2X3 A 2^max -^min 

with 

a^max = max{a;i, x^, x^} , Xmin = min{a;i, X2, x^] . (3.3) 

At the Born level, all three partons are in the event plane so that D = A3 = 0. With 
account of secondary partons /c^, radiation out of the event plane is generated and 
one gets A3 > 0. At the same time, the hard quark, antiquark and gluon momenta 
acquire recoils and move out of the event plane. In the near-to-planar region (2.3), 
both the hard parton recoil and the secondary parton momenta can be treated as 
small (soft)^. To leading order the smallest eigenvalue A3 is given by 

^3 = E % ' (3.4) 

with kix the out-of-event plane component and cui the energy of a secondary soft 
parton i. Corrections are quadratic in the "soft" parameter kx/Q ^ D <^ 1. In 
particular, due to the presence of the energy in the denominator, the contributions 
from the recoiling hard primary partons qq, g are of second order {O (-D^)) and have 
been neglected in (3.4). The quantities A1A2 and are given, to leading order, by 
(3.2) with corrections linear and quadratic in D respectively. 



^As we have shown in [7], finite rescaling of the in-plane momenta, due to hard colhnear sphttings, 
gets absorbed into the first hard correction to the emission probabihty of soft gluons, which is then 
resummed and embodied into the radiator. Bearing this in mind, all hard parton recoils and the 
secondary parton momenta ki can be effectively treated as small. 
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3.1 Resummation of soft radiation 

The starting point for tlie parton resummation of tlie distribution (2.2) is the fac- 
torisation of soft emission from the hard parton system. The distribution for 
the emission of n soft partons from the primary qq, g system can be factorised in the 
form 

{Pa, ki...kn)c^ Ml (Pa) ■ (Pa, h . . . kn) , (3.5) 

where we have used the fact that the hard parton recoils can be neglected and the 
actual qqg momenta replaced by their Born values Pa- The first factor is the squared 
first order matrix element giving the Born distribution 



^ a, = a,{Q) , (3.6) 



dxidx2 27r {1 — Xi){l — X2) ^ 

where, from now on, xi and X2 mark the quark and antiquark invariant energy 
fractions (or vice versa), and X3 = 2 — xi — X2 is that of the hard gluon. 

The second factor Sn describes the distribution for emitting n soft partons ki 
from the hard qq, g system. This distribution can also be factorised into the product 
of independent soft emissions. The structure of this factorisation depends on the 
required accuracy. Since we aim at SL accuracy, we follow the analysis of [7] in 
which Sn is factorised at two-loop order. 

In the region (2.3), from the factorised expression (3.5) at parton level we can 
write 



^{D,yc) ^ j dxidx2Q{y3-yc) |C(Qis) • • S{D, xi,X2))> 

'^iVc) - J dxidx2 Q{y3 - yc) 



(3.7) 



dxidx2 

with y3 given by the Born expression (3.2). The first factor in the curly brackets is 
the non- logarithmic coefficient C{as) — 1 + O (ctg), which is, in general, a function 
of t/3 and D, and has a finite D — > hmit. The soft factor S{D, xi,X2) accumulates 
all logarithmic dependences on D, and is given by 

S{D,X,,X2) = Y.-^ Y[^S^iPaM---kn)Q[\3-Y.^]^ (3-8) 

n I \ I / 

(with A3 = D/27A1A2). Taking the factorised structure of Sn at two loops, this 
expression for S{D, xi,X2) is accurate to SL level, see [7]. 

Given the factorised expression for S^, in order to sum up the series in (3.8) in 
suffices to factorise the theta-function constraint by using the Mellin representation. 
We obtain 

/I 
^e^"^a(z/,a;i,a;2), (3.9) 
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where 



1 r r^^k- 

a(u, ^1, ^2) = E -T / n ^ ^"'"^ • ^n{Pa: h . . . K) ^ ^-^i^,^.,^^) , (3.IO) 
n I 

The contour in (3.9) runs parallel to the imaginary axis with Re > 0. The near- 
to- planar region D corresponds to the region of large Mellin variable f ^ 1. 
We show that the radiator is given by a PT contribution and a NP correction 

7l{v, xi, X2) = T^^i^-, xi,X2) + 5Tl{v, xi,X2) ■ (3-11) 

In the next section we compute the PT contribution at SL level. In section 5 we 
compute the leading 1/Q correction including the effect of non-inclusiveness of 
the D-parameter (Milan factor). 



4. PT contribution at SL accuracy 

To obtain the PT radiator 7^^^(z/, Xi, ^2) we follow the procedure described in detail 
in [7] for the calculation of the distribution. The present case is simpler since the 
hard parton recoil can be neglected. At SL level we have 
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ri3 + - ri2 , , 



(4.1) 



where 



rah{l^-,Xi,X2) 



/ Wab{k) 


1 _ e-'A 


J TTCU 





Wab{k) 



^klab 



(4.2) 



with Wab{k) the soft distribution for the dipole ab. Here the running coupling is 
defined in the physical scheme [11] and kf^ab is the invariant transverse momentum 
of k with respect to the hard parton pair Pa, Pb, defined as 



P - 



2{Pak){kPt,) 

(PaPb) 



(4.3) 



The unity in the square bracket in (4.2) takes into account the virtual corrections. 
The contribution from the source u{k) = e^'^^'^/^'^ results from the exponentiation of 
real emissions, see (3.10). To reach SL accuracy one needs to take into account also 
the correction coming from hard coUinear parton splittings, which will be embodied 
into a redefinition of the hard scales, see later. 
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4.1 Explicit calculation of the PT radiator 

In the laboratory frame the hard parton momenta Pa and Pf, are not back-to-back 
and the expression of soft dipole distribution Wab{k) is rather cumbersome (see (4.3)). 
It is then natural to evaluate Tab starting from its expression in the centre-of-mass 
frame of the dipole Pa + Pbi 

P: = ^(1,0,0,1), P* = ^(1,0,0,-1), Ql, = 2PaP, = Q\xa+x,-l), 

(4.4) 

in which the distribution Wab{k) and the phase space are given by 

Wab{k) = , ^dK^ ^dri , 4.5 

with k'^, 7] and the squared transverse momentum, the rapidity and the azimuthal 
angle of the soft gluon, respectively. To compute the dipole radiator rab{i^) we need 
to express the source for our observable u{k) ~ e''^^'^^'^^ in the frame (4.4). While 
the variable kx is the same in the two frames, k^ — Ksm(f), the gluon energy cu in the 
laboratory frame (see (A.l) in Appendix A) is given, in terms of the dipole c.m.s. 
variables 

uj = K (^Aab cosh(r/+r/°,,) - ^A^^-l cos <j?j , (4.6) 
with Aab and the following functions of 

•^aj -^b 



Aab = -^Vx^b = \l "^"^^ - , tanh77°6 = ^° . (4.7) 

Qab \ Xa + Xb- 1 Xa + Xb 

Finally, using the frame (4.4), the rab radiator can be written in the form 
rab(u,x„X2) = J^ ——- j_^-j_Jr^(l-e = In — , (4.8) 

where r is the following function of 77 and 

sin^ 6 , ^ . 

T = ^ , ^ . (4.9) 

AabCoshiji+rPJ cos(/) 

Note that, since Qab ~ Q, to SL accuracy, we do not care for the precise upper limit 
in K as long as it is of order Q. 

The calculation of Tab is performed in Appendix B and one finds, to SL accuracy. 



Tabii^, Xi,X2) = 2 — In + 2/ — In — — - . 4.10 

Jq2/^ TT K Jq2/^2 TT \2AabQ / 
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Combining together the various pieces and recaUing the definition of Aab in (4.7) we 
find the full PT radiator 



7^P^(^/,Xl,X2) = CF[ri{u,Q^)+r2{u,Q2)] + Cat^{v,Q^) , (4.11) 
where the contribution of the hard parton #a is given by 



ra{j^,Qa)^2 — ^ In + 2/ — In — — — — 



(4.12) 



Here the hard scales Qa are 



g? = Q^ = 2(PiP2) = g'(i-a;3), 

_ 2(PiP3)(P3P2) _ ^, {l-x,){l-x,) (4.13) 

(P1P2) (1-X3) ■ 

In (4.12) we have included into the first term the rescahng factor 

6 = 6 = e-^/^ , 6 = e-^o/^^= , (4.14) 

to take into account the hard coUinear splittings of the quark and the gluon. These 
constants, as well as the precise expressions for the geometry dependent scales Qa, 
are important in order to incorporate in the final result (4.11) for the radiator all 
terms of the order of ct" In""^^ v and a" In" v. 
To first order we find 

r„(.,g„) = ^ln (4-15) 

This result corresponds, at DL level, to a half of the DL radiator in the case [7]. 
The relative factor 1/2 is due to the fact that in the case the contribution to the 
observable is linear in the angle with respect to the event plane, while it is quadratic 
in the D case. 

The first order contribution to the PT part of the distribution (3.7) is obtained 
by performing the Mellin integral in (3.9). We find 



with Ct = 20 F + Ca the sum of the colour factors — the total colour charge of the 
hard qq, g system. Here Gn is given by 

Giiivc) = J dxidx2Q{y3 - yc)-^^^^^^^^^^Gu{xi,X2) , (4.17) 



9 



with 



Gu{xi,X2) = 2Ct ln27AiA2 + 2Cf 



In 



+ In 



+ 2Ca In 



(4.18) 



We have checked that this function correctly describes the first non-logarithmic cor- 
rection by comparing with the result of the numerical program EVENT2 [12] for 
= 0.025, 0.05 and 0.1. 



5. NP calculation 

We follow the usual procedure [6] for computing the leading NP corrections, including 
two-loop order to take into account the non-inclusiveness of jet observables. We start 
from the PT expression (4.8) and then: 



we represent the coupling by the dispersive integral [5] 



Jq (k^ -|- m^)^ ' 



(5.1) 



we substitute ^/ k,"^ + m? for the momentum n in the source in (4.8); 

we take the leading part of the integrand for small k and m by linearising the 
source 



1 - e Q • 



-T , (5.2) 



since the NP part of the "effective couphng" 5aes{m) has a support only at 
small m; 

• wc multiply this expression by the Milan factor Ai, computed at two-loop 
order, to take into account effects of non-inclusiveness of jet observables. 

For the {a6}-dipole contribution this procedure gives 

(5Q;eff(m) —d f d^K + 



m 



5r.,{u,x,,x,)^uMjdm ^^^,J .(.2 + ^2) g 7^^^"^^^' 

(5.3) 

where r is the function of Xi,X2 defined in (4.9). It is related with the ratio of 
transverse momentum to energy and therefore decreases exponentially in rapidity. 
This allows us to extend to infinity the 77 integrals by setting tjm — 00. 

Here lies the main difference with the case in which the observable was uni- 
form in rapidity so that the corresponding NP radiator involved a divergent rapidity 
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integral. There one had to keep the hard parton recoil momentum which provided 
an effective cutoff to the rapidity integral and resulted in a log-enhanced NP contri- 
bution. 

We find the leading NP correction to the {a6}-radiator 

Srab{i^, xi, X2) = I' -^gab{xi, X2) , (5.4) 
where gab is the geometry dependent function 

9ab{Xl,X2)= — , (5.5) 

J-TT ^TT Aab cosh T) - COS 0-^ A^^^^ - 1 

and a^^ is the NP parameter 

a^^ ^2mJ dm^-^^ . (5.6) 
Combining all pieces we obtain the 1/Q correction to the radiator 

5n{u) = u —A{xi,x2) , ^ = Y + 923 - j^guj ■ (5.7) 

After merging the PT and NP corrections we can write the NP parameter in the 
form 

^NP ^ |^^(^^) _ «^ _ ^„ ^ (^In £ + ^ + 1^ I , a,^ a^{Q) , 

(5.8) 



where 



1 f^' 
fJ'i Jo 



The term proportional to K accounts for the mismatch between the ms and the 
physical scheme [11] and reads 

K = nJ'2-^-]-^. (5.9) 

For the analytical expression of the Milan factor Ai see [13]. To quantify the pa- 
rameter aJ^^ we recall the expression for the NP shift At in the thrust distribution 
[6], 

NP 

At = 2Cf^ . (5.10) 



11 



6. Final result 



We are now in a position to obtain the full distribution yc) in (3.7) to the stan- 
dard accuracy. First we obtain the PT expression for the soft factor S^'^{D, xi,X2) 
by performing the MeUin integral (3.9) with the PT radiator (4.11). Then, using the 
exact result of the O (a^) matrix element calculation, we compute the first correction 
of the coefficient function C{as) in (3.7). This allows us to perform the matching of 
the resummcd and the exact result to this order. Finally, we include the leading NP 
correction from (5.7). 

6.1 PT resummed distribution 

The PT contribution is given by 



2my ^(l + 7^') ' (6.1) 

n' = -DdDn^^ {D-\xuX2) , 

This evaluation of the Mellin integral is accurate to SL level, (see e.g. [7]). In the 
calculation of the logarithmic derivative of TZ the precise expression for the hard scale 
Qa in ra{i', Qa) is beyond SL accuracy as long as it is of order Q. We can use then 
Q as a common hard scale in the SL functions r', and substitute the logarithmic 
derivative of the radiator TZ' with 

where we have explicitly used the expression (4.10). 

Again, to SL accuracy, we can expand the radiator in the exponent of (6.1) and 
write 

7^P^ (A3 \ xi, X2) = 7^P^ {D-\ xi,X2) + K ■ ln(27AiA2) , (6.3) 

where O (ctg) corrections have been neglected. We conclude by giving the SL-accurate 
expression 

S^^iD, x„ X2) = e-^"^(^-^--) . . (6.4) 

Using the expression of the PT radiator given in Appendix B and TZip given in (6.2), 
this distribution can be written as 

\nS^'^{D,xi,X2) ^ Lgi{aL) + g2{aL,xi,X2) , a = q;ms(Q) , L^-\nD, 

(6.5) 
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Figure 1: The first order coefficient function for three different values of ye- 

where only the logarithmic terms in D were kept, while the finite corrections O (as) 
dropped. This precaution is necessary in order to properly set up the procedure for 
matching the resummed expression with the exact fixed order result. The DL contri- 
bution, the first term in (6.5), does not depend on xi,X2- Therefore the dependence 
on the event geometry emerges at the level of subleading SL effects. In the first order 
we recover from (6.5) the result already obtained in (4.16). 
Finally, the resummed PT distribution is given by 



S^^{D, xi,X2) . 



(6.6) 



a{yc)dxidx2 

6.2 Matching resummed and fixed-order prediction 

Here we consider the matching with the exact result of order a^, which allows us to 
compute the first order contribution to the coefficient function C{as) in (3.7). The 
exact first order result reads 

^U^tiD, yc) = l + ^ (-Ct In^ ^ - GniVc) In + c,{D, y^)^ + O {a!) , (6.7) 

where Gu{yc) is defined in (4.17) and ci{D, yc) is a function (regular at D = 0) which 
we calculate numerically by using the four-parton numerical program EVENT2 [12]. 
This gives 



C(as) = l + -^ci(A2/c) + ... 

ZTT 



(6. 



In fig. 1 we report ci{D,yc) as a function of D for three values of yc 
We define the three matched expressions: 
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• Log-R matching: 

E^I,(Ayc)^e^-(^--=)-EP^(Ayc), (6.9) 

• R-matching: 

^7jD,y^) ^ (l + ^c^{0,y,)) ■ if^iD^y^) + g (cip,ye) - Ci(0,ye)) , 

(6.10) 

• modified R-matching: 

s^It(^, Vc) ^ (i + ^ ci(^> yc)) • sP^(i^, ye) . (6.11) 

It is easy to see that each of these expressions reproduces the first order result (6.7) 
and accounts for all terms of order a" ln™D with m > 2n — 2. To obtain also the 
a" ln^"~^ D terms one needs a second order matching (O (a^) exact matrix element 
calculation) which requires five-parton generators [14], [15]. The matching procedure 
becomes more involved and will be discussed in [16]. 

6.3 Including NP corrections 

We now consider the distribution with the NP corrections included. Using the ex- 
pression (5.7) for the NP radiator, we have 

Qm ^ fdiye;^ -L^A+n^^(u)} ^p^^n xn ^ (6 12) 



where 



5D^ — ■V{xi,X2) , P(xi,X2) = 27AiA2A(xi,X2) . (6.13) 

As for other shape observables, we have embodied the leading NP correction as a 1/Q 
shift of the argument of the PT distribution. The magnitude of the shift is determined 
by the product of the universal NP parameter and the geometry dependent function 
T>{xi,X2) with A(xi,X2) given in (5.7). 

The final expression that includes the NP correction is given (for instance, in the 
Log-R matching scheme) by 

E(D,ye) = j dxidx2e{y3 - Vc) ^^y^2^^^^ eg • S^^{D', xi, X2) , (6.14) 

where D' is the shifted variable 

The two last terms in (6.15) are relevant only at large D and have been added as to 
ensure the correct normalisation of the distribution, D' — 1 at D — 1. 
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6.4 Mean value 

We can also consider the mean value of D at fixed yc defined as 

{D)^ [\dD^^^. (6.16) 
Jo 

The integral is dominated by the region of finite D where the secondary partons 
are all hard, and out of the event plane, so that no logarithmic enhancements are 
involved here. The PT contribution (D)^^ should be then obtained from a fixed 
order calculation. We have to add to it, however, a NP correction 

{D) = {Df^ + {Df . (6.17) 

The NP correction is determined by soft (small transverse momentum) partons and 
can be obtained from the present analysis, cf . [7] . It is simply related with the average 
value of the shift as follows 

„NP r J (0) 

=- y ...,.,efe-,,)^j^^^.p(..,.,). (6.18) 

with V{xi,X2) given in (6.13)^. 
7. Numerical analysis 

We report here the numerical results for some typical values of i/c- Data are not 
yet available. The results depend on the two parameters ayig{Mz) and ao{fJ>i) (with 
III = 2 GcV) which values we fix in the range determined by the 2-jet shape analysis, 
see for instance [17]. 

In fig. 2 we plot, as a function of Q, the mean value (D) given in (6.16) for 
yc = 0.05. The leading order PT contribution is obtained from EVENT2 [12]. The 
next-to-leading order PT contribution is obtained from DEBRECEN [14]. The NP 
contribution is given by (6.18). 

We see that the NP correction is large up to LEP-II energies. Actually, the 
coefficient A{xi,X2) of the NP correction is large due to the fact that one of the 
three radiating hard partons is the gluon whose contribution to the NP radiation is 
proportional to its large colour charge, see (5.7). 

In fig. 3 we plot the mean value (D) given in (6.16) for three different values 
of yc dX Q = 91.2 GeV. The mean value decreases with yc decreasing. This is 
expected since decreasing yc one includes configurations in which the hard gluon 
become close to one of the quarks, so that the independent radiation off the gluon 
gets suppressed. In figs. 4, 5 and 6 we plot the distributions E(D,|/c) for three 

^The expression (6.18) has been checked by G.P. Salam and Z. Trocsanyi using their (unpub- 
hshed) numerical program. 
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Figure 2: The mean value {D) as a function of Q for = 0.05. Here we have taken 
a^{Mz) = 0.118 and the NP parameter ao(2GeV) = 0.44. 
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Figure 3: The mean value {D) as a function of Uc for Q = 91.2 GeV for a^{Mz) and 
ao{nj) taken at the values of the previous figure. 



values of i/c for Q = 91.2 GeV. The PT distribution is given by the first order Log-R 
matched expression in (6.9). (The second order matching will be presented elsewhere 
[16].) The full distribution is given by the expression (6.14). The PT peak is shifted 
to the right by an amount of order of the NP part of the mean values, see fig. 3. The 
shift of the PT distribution is not uniform since the xi, X2-dependent shift has been 
averaged over the Born distribution. 
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Figure 5: The distribution for Q = 91.2 GeV and Vc = 0.5 for ay^{Mz) = 0.118 and the 
NP parameter ao{iJ,i) = 0.52 at m = 2 GeV. 

8. Discussion and conclusion 

The aim of the present study is the understanding of the structure of multi-jet events 
in e~^e~. To master this field is important, especially, for hadron-hadron collisions 
where the QCD events involve more than two jets. 

This is the second near-to-planar three-jet shape variable, after Tm, studied at 
an accuracy sufficient to make quantitative predictions (SL resummation, matching 
with fixed order, NP corrections). The main difference between the D-parameter 
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Figure 6: The distribution for Q = 91.2 GeV and yc = 0.1 for ay^{Mz) = 0.118 and the 
NP parameter ao{pi) = 0.52 at fij = 2 GeV. 

and is the dependence of the observable on particle rapidity. Actually, is 
uniform in rapidity, while D is exponentially dumped. This feature has important 
consequences both at PT and NP level. 

On the PT side, the hard parton recoil explicitly enters the observable in the Tm 
case. On the contrary, in the D case it gives rise to corrections beyond SL accuracy. 
This in turn implies that the PT D distribution depends only on the geometry of 
the event (the angles between jets) and not on its colour configuration. 

In what concerns NP corrections, one finds that in the D case only soft radia- 
tion at large angles contributes to the NP shift. As a consequence, the shift is not 
logarithmically enhanced, as was the case of Tm- This is the same difference that 
occurs between 1 — T,C, M^/Q"^ and B jet shapes in the two-jet case. 

We observe that the shift for the D parameter is larger than is typical for 2-jet 
observables. This is explained by the fact that in 3-jet events we have three hard 
radiating partons whose total colour charge 2Cf + Ca is significantly larger than 
2Cf, the total colour charge of a two-jet system. As a consequence, while for 2-jet 
observables higher order NP effects come into play below the peak of the distribution, 
for the D parameter they are already relevant in the proximity of the peak. This calls 
for a deeper analysis that would address higher powers in 1/Q, for example, along 
the lines of Korchemsky-Sterman approach which was recently developed for some 
2-jet observables in [18]. The comparison with experimental data (not yet available) 
would shed light on this important point. 

Near-to-planar 3-jet events provide a new method to measure the fundamental 
QCD parameter as{Mz) and to have a further test of the universality of genuine NP 
effects in jet physics. 
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A. Born kinematics 

We select the event plane in such a way that the two Born momenta with the maxi- 
mum and minimum x-component are given by 

-fmax -^(-^max) 0, 0, 2^max) ) -fmin -^(-^min) 0, 7m, tmin) • ("^■■^) 

with Q — 2E and 



so that the thrust major is 

Tm = ^ V(1 - ^i)(l - ^2)(1 - X3) . (A.3) 

■^max 

To define these massless momenta one needs that T and Tm are restricted to the 
region 



2(1 - T)^Tr -I 



<Tm< . (A.4) 



T 

In this paper we consider T, Tm restricted to this region. 
B. The PT radiator 

Here we compute the radiator component rab{T^) given in (4.8). We change integration 
variable from 77 to r with the Jacobian 

dr , , , sin^ 6 _ , 

dy]^--K{T), K(r)^ (B.l) 

In this appendix we set Aab — > A and 77°^ rio defined in (4.7). The r^-integral can 
be divided into the positive and negative rapidity regions and one has 

rai>iiy)^r^bi^)+r-,{u), (B.2) 
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where 

The hmits are for 77 = ±770 and 77 = 77M 

A — V^^ — 1 COS0 ' ^ A cosh(r7M ± ^0) — — 1 cos 
The hmit is a function of k, given, for small k,, by 

r^(«)-^7^, Q^. = Qa.e^"" ■ (B.5) 

In (B.3) we can substitute, to SL accuracy, the factor [1 — e"'^^^] by the effective 
cutoff (C.4) and we obtain 

In order to reach SL accuracy, we have to take the running coupling at two loops 
and include the non-soft part of the parton splitting function. Before discussing this 
case, it is instructive to consider the first order term in ccg = cks{Q)- 

B.l First order result 

We change again variables q = kt and obtain 

± 2as dq d(t) d.T 



where r± is the solution of the equation 



q\ 12 sin^ (f) q 



with T^{k) the function of k given in (B.4). Again in (B.7) we take as upper limit 
the hard scale Q. For small q the r integral is 

Summing the and r"^ contributions, we get, at SL accuracy, 

dq, Qab „ as , 2 / '^Qlb 



^^-^ f ^ln% ^ f^^l • (B.9) 

Finally, using (4.1) we obtain, for fixed ag, the complete radiator Tl^'^(u) in (4.11). 
The parton #a component in (4.15) is obtained from (B.9) by taking into account the 
non-soft part of the sphtting function which gives the rescaling factors in (4-14). 
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B.2 SL result 

Here we evaluate r^^ in (B.6) with the running coupling at two loops. We consider 
the two regions of the r-integration: > Q/Vk and < Q/Vk. We then have 
two contributions 

V M l/TQ UK J 

(B.IO) 

where in the second term we made the variable change k q±/ 1^- The precise 
expression of the upper limit in the first term is beyond SL accuracy as long as of 
order Q. In the second line we have neglected terms beyond SL accuracy (a^ ln"~^ u). 

Now we can perform the 0-integration observing that only terms containing one 
power of In sin^ contribute to the radiator at SL level. Therefore we make the 
substitution 



/_ 



'''^ 2 i _ 1„ ^ - ^2 , ^2 _ QQab'^'^ab 



— In sin = — In 4 =^ Q± ^ Q± 



In conclusion the radiator may be written as 

J,l k2 2t, K j^jg, 21, 2AM ' 

One finds that the dependence on q± in the integration limits cancels out, so that, 
adding r^^ and r~^, one is left with (4.10). Combining together the various pieces 
and recalling the definition of Aab in (4.7) we find the full radiator T^'^iy) in (4.11) 
with the parton #a component given by the sum of two terms 

^ 2 l'^' d^'a,{K') ^aQa _ o a,{K') ,^ fuKQa\ (B.12) 



The hard scales are given by (4.13). The SL terms coming from the parton hard 
coUinear splitting are taken into account, see [7] and [19], simply by rescaling the 
scale Qa by the constants in (4.14) in the piece r^. To first order in ctg this 
expression coincides with the result in (4.15). 

To evaluate the two components r^,r^ it is enough to consider the following 
expression for the running coupling: 

. 2x "s A /3i tts ln(l - A)\ a^(3o 
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with 

^^^ 17X;-,Xn,-3C,n, ^ ^^^^^ 

and as = as{Q) in the physical scheme [11] related to the ms by 

as = as ( 1 + + . . . j , as = ttMs(<3)' {B.15) 

with K given in (5.9). 

We can separate the pieces depending on the geometry. We write 



a-ruL + r^'sL-^^r-^] , r„^ = r£, + rf,-ln(:^) . (B.16) 



Here the geometry dependence is explicitly expressed by the scales Qa, while the 
various r-f unctions are evaluated at the common hard scale Q. The various functions 
are 

^f^^(l(Tyn\ T^OnWMfOn n/^:^M±^ -^(P) + 

+ ^ f^(p) - ^(2P) ^ 



^ _ f^^ dK^ a,{K'^) _ 2 1 ^ _ dK^a,{K'^) _ 2 1-p 



where 



(B.17) 



p^a^^lni/, L(x) = ln(l - x) . (B.18) 
47r 



Notice that the radiator contain only DL and SL terms 



n^^^^) = LF^{p) + F2{p). (B.19) 
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C. Effective cutoff 

Performing the variable change q = kt in (B.3) the radiator is proportional to 

(C.l) 

This shows that /(g) has the form 

/(g) = F{asLq) + asG{asLq) + ... , Qig = as{Q) , Lq = ln ^ , (C.2) 

where the dots do not contribute at SL level. We consider first the contribution from 
the leading piece of /(g) 

r\Q) = e-'i) F{a,L,) . (C.3) 

Jo Q ^ ' 

To SL accuracy here we can make the substitution 

(x _ e-^*) ^ ^ ^5 _ 1^ , v^e^^v. (C.4) 
For completeness we prove this well known result (see [3]). We can write 

r\Q)^ r ^ F(a,L,) + A(g-) , (C.5) 

J q Q. 



with 



(C.6) 



for e ^ 0. We can extend the last integral to infinity and get, up to e ^ corrections 
for large i/, 

Ate).F(-«A).{l(l)^Ql±f).-}. (C7) 

Our aim now is to select g in such a way that A is beyond SL accuracy. Setting 
g = Q/V we get 

A(g) = F{a,Lq - a,d,) • /(e) , L,- = - In P . (C.8) 
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where 



l-r(l + e) e' 



(C.9) 



e 



We conclude then that 



F'{x) = F{x) 



(C.IO) 



which is neghgible within SL accuracy. We conclude then, to SL accuracy, 



For the second term asG{asLq) in (C.2) the analysis is simpler. Since it gives only 
a SL contribution the lower scale can be taken at any value of the order Q/v. We 
then conclude that to compute the radiator to SL accuracy we can use the cutoff 
substitution (C.4). 
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